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Abstract 

A  mathematical  model  for  analysis  and  prediction  of  the  performance  of  the  aluminum-air  cell  has  been  developed.  The  model  takes 
into  account  the  kinetics  of  the  anode,  cathode,  and  parasitic  reactions.  Ohmic  losses  in  the  electrolyte  and  mass  transfer  are  also  taken 
into  account.  The  model  prediction  of  cell  performance  shows  good  agreement  with  experimental  data.  The  mathematical  model  provides 
detailed  information  about  cell  performance  for  a  wide  range  of  operating  and  design  parameters.  For  better  cell  performance,  our  model 
studies  suggest  the  use  of  higher  electrolyte  flow  rates,  smaller  cell  gaps,  higher  conductivities,  lower  parasitic  current  densities  and 
operation  at  moderate  current  density.  From  our  analysis,  we  have  determined  that,  in  an  aluminum  air  cell,  only  the  activation  and  ohmic 
overpotential  are  important. 
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1.  Introduction 

An  aluminum-air  battery  system  has  the  potential  to  be 
used  to  produce  power  to  operate  cars  and  other  vehicles 
[1,2].  There  are  many  parameters  involved  in  the  design  and 
operation  of  the  aluminum-air  battery  system.  To  optimize 
design  and  operation  of  an  aluminum-air  battery  system,  it 
is  essential  to  develop  a  mathematical  model  that  can  predict 
the  aluminum-air  cell  performance.  In  our  previous  paper, 
we  briefly  provided  the  cell  performance  model  equations 
[3].  In  this  paper,  a  mathematical  model  for  analysis  and 
prediction  of  the  performance  of  the  unit  cell  has  been  de¬ 
veloped  in  detail.  The  effects  of  mass  transfer,  migration, 
gas  evolution,  kinetics  of  electrodes,  parasitic  reaction  and 
ohmic  losses  in  the  electrolyte  are  considered.  By  intro¬ 
ducing  selectivity,  the  algorithm  for  modeling  calculations 
is  simpler.  MATLAB  software  programs  were  used  for 
these  calculations.  This  model  will  be  used  to  predict  cell 
performance  (current  density-cell  voltage  curve),  the  ef¬ 
fect  of  electrolyte  flow  rate  on  current  density  distribution, 
the  effect  of  cell  gap  on  current  density,  and  the  effect  of 
electrolyte  conductivity  on  current  density.  It  also  is  used 
to  predict  the  effect  of  parasitic  current  on  selectivity,  ac¬ 
tivation  overpotential  and  ohmic  loss  at  different  cell  gaps, 
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gas  fraction  variation  with  cell  height,  and  the  effect  of  gas 
fraction  on  electrolyte  conductivity.  The  model  also  predicts 
the  anode  surface  concentrations  of  OH-  and  Al(OH)4_, 
and  mass  transfer  boundary  layer  thickness. 


2.  Problem  definition,  geometry  and  assumptions 

The  cell  is  made  up  of  two  plane  electrodes  with  cell  gap 
(S)  and  height  (//),  where  the  height  is  much  greater  than 
the  cell  gap  (Fig.  1  shows  the  geometry).  The  electrolyte 
enters  from  the  bottom  in  laminar  flow  with  developing 
boundary  layer  (i.e.  the  boundary  layer  thickness  increases 
with  height).  The  electrochemical  reactions  occur  only  on 
the  electrode  surfaces.  No  crystallization  reaction  occurs  in 
the  cell. 

At  the  anode,  the  main  reaction  is 
A1  +  40H"  -+  Al(OH)4"  +  3e“. 

The  parasitic  (undesired)  reaction  at  the  anode  is 
A1  +  3H20  +  OH"  |H2  +  Al(OH)4". 

At  the  cathode,  the  reaction  is 
02  +  2H20  +  4e“  -+  40H". 
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Nomenclature 


A,  b 

Ci 

CiO 

cib 

ci, ref 

^bubble 

Di 

Dr 

E 

Eoi 

Eeq 

F 

fly ) 

G 

H 

h 

Gv 

ic 

finO 

h 

ip() 


N 

Ni 

P 

R 

Rf 

Ri 

Re  y 


S 

Sm/p 

Sc 

t 


empirical  cathode  kinetic  constants 
concentration  of  species  i  (mol/cm3) 
surface  concentration  of  species  i 
(mol/cm3) 

bulk  concentration  of  species  i 
(mol/cm3) 

reference  concentration  of  species  i, 
in  this  article,  c;,ref=Qb  (mol/cm3) 
diameter  of  bubbles  (cm) 
diffusion  coefficient  of  species  i 
(cm2/s) 

diffusion  coefficient  of  the  limiting 
reactant  (cm2/s) 

empirical  cathode  kinetic  constants 
equilibrium  potential  of  electrode  1  V 
cell  equilibrium  voltage  (V) 

Faraday’s  constant  (96,500 C/eq.) 
gas  fraction  as  a  function  of  height 
gravitational  acceleration  (cm/s2) 
electrode  height  (cm) 
anode  current  density  (mA/cm2) 
average  current  density  (mA/cm2) 
cathode  current  density  (mA/cm2) 
main  anode  reaction  current 
density  (mA/cm2) 
exchange  current  density  of  the 
main  anode  reaction  (mA/cm2) 
parasitic  current  density  (mA/cm2) 
exchange  current  density  of  the  anode 
parasitic  reaction  (mA/cm2) 
electron  transfer  number  in 
anode  main  reaction 
electron  transfer  number  in 
anode  parasitic  reaction 
electron  number  per  mol 
hydrogen  molecule 
flux  of  species  i  (mol/(cm2  s)) 
pressure  (dyn/cm2) 
ideal  gas  constant  (J/(mol  K)) 
empirical  cathode  kinetic  constant 
reaction  rate  of  species  i  (mol/(cm3  s)) 
Reynolds  number  as  a  function  of  y 
stoichiometric  coefficient  of  anode 
main  reaction 

stoichiometric  coefficient  of  anode 
parasitic  reaction 
cell  gap  (cm) 

selectivity  of  main  reaction  over 
parasitic  reaction 
Schmidt  number  of  species  i 
in  the  electrolyte 
time  (s) 


T 

temperature  (K) 

ut 

mobility  of  species  i  (cm2  mol/(J  s)) 

< 

bubble  rising  velocity  (cm/s) 

V 

fluid  velocity  (cm/s) 

Fcell 

cell  voltage  (V) 

x,  y 

axial  position  (cm) 

X 

cathode  extension 

Zi 

charge  valence  of  ions  i 

Greek  letters 

^ m  5  Yim 

empirical  anode  kinetic  constants 

Yip 

empirical  parasitic  kinetic  constants 

Si 

diffusion  layer  thickness  of  species  i  (cm) 

$ic 

diffusion  layer  thickness  of  species  i 
due  to  convection  (cm) 

Sib 

diffusion  layer  thickness  of  species  i 
due  to  bubble  effect  (cm) 

0 

solution  potential 

01,  02 

electrode  1  (anode)  potential  and 
electrode  2  (cathode)  potential 

j  ^b 

solution  potential  just  outside  the  diffusion 
layer  at  electrode  1  (anode)  and  solution 
potential  just  outside  the  diffusion  layer  at 
electrode  2  (cathode) 

m>  m 

overpotenial  of  electrode  1  (anode) 
and  2  (cathode) 

?7al 

activation  overpotential  of  electrode  1 
(anode)  (V) 

Ocl 

concentration  over  potential  of  electrode  1 
(anode)  (V) 

0  m 

activation  overpotential  of  anode  (V) 

»?p 

activation  overpotential  of  parasitic 
reaction  (V) 

K 

conductivity  of  electrolyte  at  gas 
fraction /(y)  (S/cm) 

K0 

conductivity  of  electrolyte  at  zero  gas 
fraction  (S/cm) 

l1 

viscosity  of  electrolyte  (g/(cm  s)) 

P 

density  of  electrolyte  (g/cm3) 

Pgas 

density  of  hydrogen  gas  (g/cm3) 

3.  Fundamental  equations 

3.1.  Transport  in  electrolyte  solutions 

The  laws  of  transport  in  dilute  electrolyte  solutions  have 
been  known  for  many  years  and  have  been  discussed  in 
detail  elsewhere  [4].  The  flux  of  a  species  due  to  migration 
in  an  electric  field,  diffusion  in  concentration  gradient,  and 
convection  with  the  fluid  velocity  is: 


Ni  =  -ZiUiFci\/0  -  DiVci  +  vci 


(1) 
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Fig.  1.  Schematic  of  the  geometry  of  an  aluminum-air  cell:  <$*,  the 
boundary  layer  thickness  of  species  /;  H,  cell  height;  S,  cell  gap;  H  5; 
v,  electrolyte  velocity;  v  is  measured  from  the  electrode  surface  and  y  is 
measured  from  the  cell  entrance. 


A  material  balance  for  a  small  volume  element  in  the 
electrolyte  is: 

5  =  +  Ri  (2) 

ot 

In  our  model,  there  are  no  reactions  in  the  electrolyte 

(Ri  =  0). 

To  a  very  good  approximation,  the  solution  is  electrically 
neutral, 

J2  ZiCi  =  o,  (3) 

i 

except  in  the  diffuse  part  of  the  double  layer  very  close  to 
an  interface.  The  current  density  in  the  electrolyte  solution 
is  due  to  the  motion  of  charged  species: 

i  =  FJ2ziNi.  (4) 

i 

These  equations  provide  the  basis  for  the  analysis  of 
transport  in  the  electrolyte.  The  flux  relation  (Eq.  (1)) 
defines  transport  coefficients — the  mobility  u\  and  the  dif¬ 
fusion  coefficient  Dz  of  an  ion  in  a  dilute  solution.  For  flow 
of  the  electrolyte,  the  fluid  velocity  is  determined  from  the 
Navier-Stokes  equation: 

p  =  S/p  -  jlS72v  +  pg.  (5) 

Eqs.  (1)— (5)  can  be  combined  to  determine  the  concen¬ 
tration  of  each  species: 


where 

K  =  F2^iUiCi  (8) 

i 

is  the  solution  conductivity.  Multiplication  of  Eq.  (6)  by  ZiF 
and  addition  over  i  gives  an  equation  for  conservation  of 
charge: 

— V  •  i  =  0  =  kW20  +  (V/c)  •  V0  +  'ZjDj^Cj.  (9) 

i 

The  Nernst-Einstein  relation  relates  ionic  diffusion  coeffi¬ 
cients  and  mobilities,  at  least  approximately: 

A  =  RTui .  (10) 


We  take  all  the  diffusion  coefficients  to  be  of  roughly  the 
same  magnitude,  but  all  of  them  to  be  small  in  the  sense  that 
the  Peclet  number  Pe  is  large: 


UL 

Pe  =  - 

Dr 


(11) 


where  U  is  the  characteristic  velocity,  L  the  characteris¬ 
tic  length,  and  Dr  is  the  characteristic  diffusion  coefficient 
(usually  taken  to  be  that  of  the  limiting  reactant).  Note  that 
when  Dr  goes  to  zero  the  Pe  becomes  infinite. 


3.2.  Electrolyte  bulk  medium 


The  region  of  the  electrolyte  outside  the  boundary  layer  is 
called  the  bulk  region.  It  is  widely  known  that  when  the  Pe  is 
large,  mass  transfer  by  convection  dominates  over  diffusion 
except  in  a  thin  region,  known  as  the  diffusion  layer,  near 
the  reaction  surface  (the  electrode  in  this  case).  For  large  Pe 
(in  the  bulk  medium),  Eq.  (6)  reduces  to 

9c  ■ 

-1  +  v  ■  Va  =  0.  (12) 

ot 

That  is,  the  concentration  of  a  fluid  element  is  constant  as 
it  moves  through  the  solution.  In  most  cases  the  appropriate 
solution  to  Eq.  (12)  in  the  bulk  medium  is 

ci  =  c/b  (13) 

and  all  concentrations  have  their  bulk  values. 

For  the  region  outside  the  diffusion  layer,  Eq.  (13)  ex¬ 
presses  the  solution  for  the  concentrations.  It  is  still  neces¬ 
sary  to  solve  for  the  potential  by  means  of  Eq.  (9),  which 
in  the  bulk  solution  reduces  to 


V  •  =  0. 


(14) 


9c  ’ 

—2  +  i iVcj  =  DiV2Ci  +  ZiUiFciS72<P  +  ZiUiF(Vci)  ■  S/0, 
ot 

(6) 

The  current  density  can  be  expressed  as 
i  =  -kV0  -  Fj^ZiDiVa,  (7) 


3.3.  Electrolyte  diffusion  layer 

Diffusion  cannot  be  neglected  in  the  layer  near  the  reac¬ 
tion  surface.  But  other  simplifications  are  still  possible.  On 
account  of  the  thinness  of  the  diffusion  region,  effects  of 
curvature  can  be  neglected.  We  adopt  the  boundary  layer 
coordinates  where  v  is  the  normal  distance  from  the  surface 
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and  y  is  measured  along  the  surface  from  the  entrance.  In 
the  diffusion  layer,  Eq.  (6)  simplifies  to 


the  boundary  conditions  for  the  concentrations  in  the  diffu¬ 
sion  layer  are 


dot 

~dt 


+  v 


dCi 

dx 


+  Vy 


dCi 

dy 


dx2 


+  ZiUjF 


d20  3 C[  30 
1  dx2  dx  dx 


(15) 


We  can  ignore  the  derivatives  with  respect  to  y  compared 
to  the  derivatives  with  respect  to  x.  One  more  simplification 
is  possible.  We  assume  that  the  Schmidt  number  Sc  is 

&=-J-»l,  (16) 

14  R 

where  v  is  the  kinematical  viscosity  of  the  fluid.  The  Sc 
is  on  the  order  of  1000  for  the  electrolytic  systems  of  in¬ 
terest  here.  With  the  assumption  that  the  Sc  is  large,  the 
diffusion  layer  is  thin  even  when  compared  with  any  hydro- 
dynamic  boundary  layer,  which  may  be  present,  and  within 
a  two-dimensional  diffusion  layer  the  velocity  components 
can  be  represented  as 

vx  =  -\x20(y)  and  Vy=xfi(y),  (17) 


where  fi(y)  is  the  velocity  derivative  at  the  solid  wall, 
/3  =  dvy/dx  at  x  =  0,  and  the  prime  denotes  the  derivative 
with  respect  to  y.  These  are  the  first  terms  in  expansions  of 
the  velocity  in  x  and  satisfy  the  continuity  equation: 


dVx 

dx 


+ 


dVy 

dy 


(18) 


With  this  approximation,  the  diffusion  layer  equation  for 
the  concentrations  is 


dci 

~dt 


dCi 

dx 


+  x/3 


dci 

dy 


d2Cj 

dx 2 


+  zmF 


dz0  dci  d0 
1  dx2  dx  dx 


(19) 


These  equations  (one  for  each  species,  OH-,  Al(OH)4-) 
are  to  be  solved  along  with  the  equation  of  electro-neutrality 
(Eq.  (3))  and  certain  boundary  conditions  that  are  yet  to  be 
discussed. 


Ci  ^  Cibi  % 

->  5m 

(21) 

,T  Sjmlm 

Vm  = 

nmF 

at  x  =  0, 

(22) 

The  parasitic  reaction  at  the  anode  is  expressed  as 

E  s/p M/  o.  (23) 

i 


We  obtain  the  stoichiometric  coefficient  from  the  reaction 
Eq.  (23)  which  represents  the  net  parasitic  corrosion  reac¬ 
tion  and  we  obtain  the  electron  transfer  number  rip  from 
the  oxidation  corrosion  reaction.  Then  the  parasitic  flux  of 
species  i  is 


siplp 

ripF 


(24) 


So,  the  total  flux  of  species  i  is 


IS iy  —  W'm  +  IS iD  — 


siplp 


ip 


at  x  =  0. 


nmF  cipF 
If  we  define  selectivity  of  im  over  ip  as  5m/p, 


(25) 


5m  /  p  — 


'm 


(26) 


Then  Eq.  (25)  can  be  changed  into 


IS iy  —  N[  m  T  IS  jj)  — 


Si 


i  m 


sip 


ip 


nmF  ?ipFSm/p 


m 


at  x  =  0. 
(27) 


4.2.  Boundary  conditions  for  potentials  in  the  bulk  medium 

The  potentials  on  the  electrode  are  shown  in  Fig.  2.  0\ 
is  the  potential  at  the  anode  and  02  is  the  potential  at  the 
cathode.  The  cell  voltage  is  the  difference  between  these  two 
potentials 

Vcell  =  02  —  0 1,  (28) 


4.  Boundary  conditions 

4.1.  Flux  and  concentration  on  the  boundaries 

The  differential  equations  describing  the  electrolyte 
solution  require  boundary  conditions.  The  most  complex  of 
theses  concerns  the  kinetics  of  electrode  reactions.  A  single 
electrode  reaction  can  be  written  in  symbolic  form  as 

^  ^  Sim^Slf'1  — >  >  (20) 

i 

where  Sim  is  the  stoichiometric  coefficient  of  species  i  and 
M i  is  a  symbol  for  the  chemical  formula  of  species  i.  Then 


e 
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Fig.  3.  Visualization  of  a  simplified  model  for  anodic  oxidation.  <£a  is 
the  boundary  condition  for  the  potential  equation:  V  •  k57<3>  =  0  in  the 
bulk  medium  (adapted  from  [5]). 


Combining  Eqs.  (28)  and  (29)  with  Eqs.  (31)  and  (32) 
gives  the  potential  at  the  anode: 

^a  =  ^eq  Ece]]  T]\  (33) 

and  the  potential  at  the  cathode: 

d>c  =  ->/2-  (34) 


5.  Additional  conditions 


5.1.  The  concentration  and  activation  overpotentials 


For  a  solution  with  an  excess  of  supporting  electrolyte, 
the  concentration  overpotential  in  the  electrolyte  diffusion 
layer  is  [4] 


=  E 

i 


Si  Dj_ 

nFui 


In 


(35) 


The  surface  or  activation  overpotential  as  defined  above 
should  depend  only  on  the  reaction  rate  and  the  concentra¬ 
tions  at  the  electrode  surface: 


The  equilibrium  cell  voltage  (Nernst  voltage)  is 

q  =  E02  ~  E o\.  (29) 

The  boundary  conditions  are  affected  by  the  sum  of  both 
the  activation  and  the  concentration  overpotential.  At  the 
anode: 


7?a  =  fUx  (x  =  0),  c/o,  •  •  •  ].  (36) 

This  can  be  obtained  by  the  modified  Butler- Volmer  or 
modified  Tafel  equations  (see  Section  5.4). 

5.2.  Concentration  profile 


0  =  OaUx  (x  =  0),  C/o,  . . .  ]  +  rjcUx  (x  =  0),  C/o,  C/,  .  .  .  ]. 

(30) 

This  equation  can  also  be  written  for  the  cathode  at 
x  =  S. 

The  boundary  condition  depends  on  the  local  current  den¬ 
sity  flowing  through  the  diffusion  layer  and  the  double  layer. 
Fig.  3  visualizes  a  simplified  model  for  anodic  oxidation 
used  to  develop  the  boundary  conditions  for  the  potential 
Eq.  (14). 

The  potential  just  outside  the  diffusion  layer  at  the  anode 
side  is 

—  Eoi  rj  i,  (31) 

where  Eoi  is  equilibrium  potential  of  the  anode,  iji  is  the 
overpotential  just  outside  the  diffusion  layer.  The  boundary 
layer  includes  both  double  layer  and  the  diffusion  layer.  In 
an  electrochemical  system,  there  at  least  two  electrodes,  the 
boundary  conditions  for  the  bulk  solutions  at  the  cathode 
side  is 

@c  =  @2  ~  E02  ~  02’  (32) 

Written  under  this  form  the  boundary  conditions  are  not 
practical  for  use  because  only  voltage  differences  are  of 
importance. 


The  concentration  profile  is  linearized  and  the  whole 
problem  is  divided  into  two  parts: 

•  a  thin  layer,  called  the  stagnant  diffusion  layer  with  a 
thickness  <$/; 

•  the  bulk  of  the  solution  where  no  concentration  gradients 
exist  (see  Fig.  4). 

5.3.  Surface  concentration  on  the  anode 

In  the  diffusion  layer,  the  flux  of  OH-  and  Al(OH)4~ 
due  to  the  migration  and  diffusion,  can  be  simplified  to  one 
dimension  in  the  y  direction  for  each  component.  Assuming 


Fig.  4.  The  concentration  profile  near  the  electrode.  The  Nernst 
diffusion-layer  is  obtained  by  extrapolating  the  linear  part  of  the  concen¬ 
tration  change  to  the  concentration  of  the  bulk  solution. 
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a  one  dimension  linear  concentration  profile  and  Ohm’ s  law 
holds: 


(Cib  GiO))x 
Ci  =  CiQ  H - - - 


and 


-V0  = 


bn 

*()’ 


Eq.  (1)  becomes: 


Ni  =  ZiUiFci 


bn 

KO 


Ji 


(37) 


and  let  q  =  (qb  +  qo)/2,  and  combining  Eq.  (37)  with 
Eq.  (27),  yields  the  concentration  at  the  surface  as  a  function 
of  y  as 


CjQ  =  1  +  KOm//c0)(E/27?7)] 

Qb  1  -  [^(/m//c0)(E/2/?7)] 

[fcm/%f)  T  (^zp/^p^£m/p)]hn 
”  Di(cib/8i)[l  -  8i(im/ico)(F/2RT)Y  (  } 

where  <$z  is  the  diffusion  layer  thickness  for  species  i  and  is 
estimated  by  the  convective  electrolyte  flow  and  the  micro 
convection  of  gas  evolution  (see  Section  5.5).  is  defined  as 
the  distance  from  the  electrode  where  q  =  qb,  by  assuming 
a  linear  concentration  profile  of  species  i. 

Eq.  (38)  can  be  written  in  terms  of  the  corrosion  current 
density,  /p,  and  the  selectivity  Sm/p: 


CiQ  1  +  \8i(ipSm/p/ kq)(F/2RT)] 

Qb  =  1  -  mvSm/v/K0)(F/2RT)] 

[fcm^m/p/^mf)  T  fe'p/^pE)]ip 
~~  Di(cib/8i)[l  -  8i(ivSm/v/K0)(F/2RT )] '  (  } 

Our  equations  to  predict  surface  concentrations  of  OH- 
and  Al(OH)4-  are  different  from  Chan  and  Savinell  [6].  The 
difference  is  in  the  coefficients  of  these  equations.  We  use 
Eq.  (20)  (main  reaction)  and  Eq.  (23)  to  accurately  determine 
these  coefficients.  Si/n  is  the  stoichiometric  coefficient  over 
the  electronic  transfer  number. 

The  main  reaction  in  Eq.  (20)  form  is 


A1  +  40H-  -  Al(OH)4-  =  3e“. 


Therefore,  Sim/nm  =  4/3  for  OH  and  Sim/nm  =  — 1/3  for 
Al(OH)4-. 

The  parasitic  reaction  in  Eq.  (23)  form  is 

A1  +  3H20  +  OH"  -  |H2  -  Al(OH)4"  -*  0. 

Therefore,  qp/np  =  1/3  for  OH-  and  q-p/np  =  — 1/3  for 
Al(OH)4~  where  qp  comes  from  the  net  parasitic  (corro¬ 
sion)  reaction  and  nv  comes  from  the  oxidation  corrosion 
reaction. 


5.4.  Electrode  kinetics 

5.4.1.  Main  anode  reaction  kinetics 

Polarization  data  for  aluminum  anodes  have  been  reported 
by  Rudd  [7].  A  simple  modified  Butler- Volmer  form  was 
used  in  this  modeling  work.  The  data  of  Rudd  was  fitted  to 
parameters  of  the  following  expression: 


where  r)m  =  —  r]ci  ~  Eq\.  The  subscript  m  stands 

for  the  main  reaction. 


5.4.2.  Parasitic  anode  reaction  kinetics 

The  parasitic  reaction  on  aluminum  can  be  written  as 


h — *p°n 


CjQ 
ci,  ref 


Ypi 


exp 


ap  F 
RT 


(41) 


where 


Op  —  ?7la  Eq i  Eop. 


(42) 


5.4.3.  Cathode  kinetics 

For  the  porous  air  cathode,  a  modified  Tafel  equation  was 
used: 

0c  =  a  +  Mog(|ic|)  +  e\ic\.  (43) 


5.5.  Diffusion  layer  thickness 


The  diffusion  layer  thickness  is  estimated  by  taking  into 
account  convection  from  electrolyte  flow  and  microconvec¬ 
tion  effects  from  gas  evolution  (see  also  [6]). 

The  position  dependent  boundary  layer  thickness  due  to 
convection  effects  is  expressed  as  [8] 

Sic  =  1.08  Sc-l/3Re-l/2y°-95S0-05.  (44) 


In  a  stagnant  fluid,  the  diffusion  layer  thickness  can  be 
correlated  to  the  bubble  generation  rate  by  the  equation  [9] : 

Sib  =  0.0 103/“°-36.  (45) 


The  combined  effect  of  convection  and  microconvection 
on  the  diffusion  layer  thickness  can  be  estimated  as  [10] 


Si  = 


8ir8 


ic°ib 


(4 +  4) 1/2 


zb 


(46) 


A  calculated  value  of  <$;( OH-)  at  a  cell  voltage  of  1.3  V, 
height  of  17  cm  and  a  cell  gap  of  0.2  cm  is  about  6.5  x 
10-3  cm,  3.2%  of  the  cell  gap. 


5.6.  Electrolyte  conductivity 


The  presence  of  the  gas  bubbles  in  the  cell  gap  decreases 
the  electrolyte  conductivity,  and  this  effect  can  be  estimated 
by  the  Bruggemann  equation  (see  also  [6]): 

k(v)  =  k()[1  -  fly)]'  5  ,  (47) 
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wher e/(y)  is  the  axial  position  y  dependent  gas  fraction  of  the 
electrolyte,  fly)  can  be  estimated  by  the  following  equation 
[11]: 


^gas 

HnFv 


(48) 


where  v'b  is  bubble  rise  velocity,  n  is  electrons  trans¬ 
ferred/mol  gas,  and  yOgas  is  gas  density. 


5.7.  Bubble  rise  velocity 


In  a  stagnant  solution,  after  formation,  a  bubble  rapidly 
accelerates  to  its  terminal  velocity,  v b.  The  value  of  r>b  is 
determined  by  the  balance  between  the  buoyant  raising  force 
and  the  drag  force.  While  it  is  easy  to  calculate  the  buoyancy 
force  for  a  bubble,  the  drag  force  varies  with  bubble  size.  For 
small,  spherical  bubbles,  the  drag  force  can  be  calculated, 
and  when  combined  with  the  buoyancy  force,  yields  Stoke’ s 
law: 


..  PSd bubble 

Vh  =  - 

18/x 


(49) 


where  p  is  the  density  of  the  fluid,  ^bubble  the  diameter  of 
the  bubbles,  and  /x  is  the  viscosity  of  the  fluid. 

If  the  velocity  of  the  fluid  is  v,  the  rise  velocity  of  the 
small  bubble  will  be: 


v'b  =  Vb  +  V 


(50) 


6.  Cell  performance  model 

6.1.  Cell  performance  model:  summary 

Eqs.  (26),  (33),  (34),  (38)  and  (40)-(46)  form  the  set  of 
equations  for  our  cell  performance  model. 

The  model  includes  the  effects  of  mass  transfer,  migra¬ 
tion,  gas  evolution,  and  the  kinetics  of  the  electrodes.  The 
mathematical  model  provides  detailed  information  about 
cell  performance  for  a  wide  range  of  operating  and  design 
parameters.  In  summary,  the  electrolyte  has  been  broken 
into  two  regions:  bulk  medium  and  diffusion  layer.  If  we 
combine  the  diffusion  layer  and  double  layer  as  one  bound¬ 
ary  layer,  the  cell  performance  of  the  aluminum-air  cell  can 
be  modeled  by  the  potential  equation: 

V  •  a;V0  =  0. 

boundary  conditions: 

=  ^eq  Vcell  ~  0l-> 

=  -rj  2, 

and  additional  Eqs.  (38)  and  (40)-(46)  which  are  used  to  cal¬ 
culate  r]\  and  t]2  (see  Section  5).  The  calculation  procedures 
are  outlined  in  Section  6.3.  We  determine  current  density  at 
position  v  and  average  current  density  of  the  cell  at  a  fixed 
cell  voltage. 


6.2.  Model  assumptions 

1.  Two  electrodes  are  in  parallel. 

2.  The  cell  is  sufficiently  long  (H  S)  so  that  entrance 
effects  on  the  cell  performance  curve  are  neglected. 

3.  Simple  Ohm’s  law  is  used  in  the  bulk  layer  to  predict  the 
potential:  —  V<£>  =  (<£>a  —  d>c)/S  =  im/ k  (one  dimension). 

With  these  assumptions,  the  cell  performance  modeling  is 
turned  into  a  one-dimensional  problem. 

6.3.  Model  computational  algorithm 

For  a  fixed  cell  voltage,  the  following  procedures  are  used 
to  calculate  current  density  at  position  v  and  average  current 
density  of  the  cell: 

1.  Assuming  a  current  density  /m,  c/o  =  Qb  (no  concentra¬ 
tion  polarization)  and  k  =  ko  (no  gas  effect). 

2.  Calculate  rial  by  Eq.  (40),  rjal  by  Eq.  (43),  ijp  by  Eq.  (42) 
and  ip  by  Eq.  (41). 

3.  Calculate  8i  by  Eqs.  (44)-(46),  fly)  by  Eq.  (48),  /c(y)  by 
Eq.  (47). 

4.  Calculate  c;o/Cb  by  Eq.  (38)  and  by  estimating  V0(= 
(0a  —  <pc)/ s  —  lm/ k). 

5.  Recalculate  rial  by  Eq.  (40),  0a2  by  Eq.  (43),  Op  by 
Eq.  (42)  and  ip  by  Eq.  (41),  and  calculate  5m/p  by 
Eq.  (26). 

6.  Calculate  voltage  of  the  cell,  V. 

Calculate  V  by  Eqs.  (33)  and  (34)  and  (<£a  —  <PC)/S  =  im/K 
(Ohm’s  law): 

Case  1:  V  >  VCq\\,  increase  im,  go  back  to  step  2. 

Case  2:  V  <  VCq\\,  decrease  zm,  go  back  to  step  2. 

Case  3:  V  =  Vce\\,  the  resulted  im  is  the  current  density 
at  position  y. 

7.  Calculate  average  current  density 

Efiyii n 

The  parameters  used  to  calculate  cell  performance  are 
tabulated  in  Table  1 . 

If  we  determine  the  cathode  parameters  from  the  Yardney 
(AC  series)  cathode  [12]  for  oxygen  reduction,  the  coeffi¬ 
cients  in  Eq.  (43)  listed  in  Table  2. 


7.  Model  calculations 

7.1.  Model  verification 

Fig.  5  shows  the  model  predicted  cell  performance  curve 
(curve  1)  compared  to  experimental  data.  The  solid  curves 
is  plotted  using  the  Algorithm  in  Section  6.3.  The  circles 
are  experimental  data  at  temperature  of  333  K  (see  [6])  with 
electrolyte  velocity  of  1.7  cm/s  and  a  cell  gap  of  0.14  cm.  In 
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Table  1 

Parameters  used  in  the  cell  performance  modeling  [6] 
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Parameters 

Anode  kinetic  parameters 

Main 

im0  =  13.71  mA/cm2,  am  =  0.07956,  ym  =  0.5,  Eoi  =  —2.4403  V  (from  Eq.  (40)) 

Parasitic 

ip0  =  11  mA/cm2,  ap  =  —0.0591,  yp  =  1,  Eop  =  0.9058  V  (from  Eqs.  (41)  and  (42)) 

Cathode  kinetic  parameters 

See  Table  2  (from  Eq.  (43)) 

Substance  physical  properties 

c0H_  b  =  5  x  10-3  mol/cm3,  D0 H-  b  =  5.26  x  10-5  cm2/s,  cA10H-  b  =  5  x  10~4  mol/cm3, 
DAioh~  b  =  10~5  cm2/s>  P  =  1.15  g/cm3,  /x  =  0.008  g/(cms),  pgas  =  7.32  x  10~8  g/cm3, 
^bubble a  =  2.6  x  10-3  cm,  kq  =  0.8  s/cm 

Operating  conditions 

Temperature  =  333  K 

Constants 

F  =  96,500  C/eq.,  R  =  8.314  J/(molK),  g  =  980  cm/s2 

a  From  [9] . 

Table  2 


Cathode  kinetic  parameters  for  the  modified  Tafel  equation  (Eq.  (43)) 


Cathode 

Kinetic  parameters 

a 

b 

e 

AC65  Yardney  (Ag) 

-0.2856 

-0.0316336 

-6.631  x 

10“4 

AC75  Yardney  (CoTMPP) 

-0.1945 

-0.02901276 

-7.402  x 

10“4 

AC78  Yardney  (Pt) 

-0.2118 

-0.0255886 

-6.706  x 

10“4 

Chan  and  Savinell 

-0.29289 

-0.025096 

-2.03429 

x  10“4 

curve  1,  the  cathode  kinetic  parameters  are  from  Chan  and 
Savinell.  The  agreement  between  the  model  and  the  data  is 
good  at  current  densities  from  180  to  500mA/cm2. 

In  curve  2,  the  cathode  kinetic  parameters  are  from  Yard- 
ney  AC78.  From  the  kinetic  parameters  of  Yardney  cathode 
AC78,  we  can  see  that  the  ohmic  resistance  term  e  (absolute 
value)  is  bigger  than  the  one  given  by  Chan  and  Savinell, 


so  at  higher  current  density,  the  Yardney  AC78  cathode  has 
higher  polarization  thus  the  cell  voltage  is  lower  at  the  same 
current  density. 

7.2.  Effect  of  electrolyte  flow  rate  on  current  density 
distribution 

Fig.  6  shows  the  effect  of  electrolyte  flow  rate  (velocity) 
on  current  density  distribution  along  the  cell.  Because  of  the 
evolution  of  hydrogen  gas  in  the  electrolyte,  the  ohmic  re¬ 
sistance  will  increase  along  the  cell  (resulting  in  a  decrease 
of  current  density  along  the  cell)  if  the  electrolyte  flow  rate 
is  not  high  enough  to  remove  the  bubbles.  With  increase 
in  the  electrolyte  flow  rate,  the  current  distributions  along 
the  cell  become  closer  to  uniform,  which  suggests  that 
operating  at  a  high  flow  rate  (8  cm/s),  the  gas  effect  will  be 
alleviated. 


Fig.  5.  Cell  performance  of  an  aluminum  air  cell  with  pure  aluminum  anode.  Curve  1:  the  cathode  kinetic  parameters  are  from  Chan  and  Savinell; 
2:  the  cathode  kinetic  parameters  are  from  Yardney  AC78. 


curve 
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Fig.  6.  Effect  of  flow  rate  on  current  distribution  along  the  cell.  Cell  gap  =  0.2  cm. 


Table  3 


Velocity  and  Reynolds  number  relationship 


Velocity  (cm/s) 

Re 

1 

56 

3 

167 

5 

279 

8 

446 

14 

780 

18 

1003 

Cell  gap  —  0.2  cm. 


The  Reynolds  number  (Re)  increases  with  flow  rate. 
The  relationship  of  Re  and  velocity  is  tabulated  in 
Table  3.  The  Re  number  indicates  laminar  flow  over  a 
wide  range  of  velocities  because  of  the  narrow  channel 
used. 


7.3.  Effect  of  cell  gap  on  current  density 

The  effect  of  cell  gap  on  the  current  density  at  different 
constant  voltages  is  shown  in  Fig.  7.  To  alleviate  the  effect 
of  gas  bubbles  in  the  electrolyte,  a  flow  rate  of  8  cm/s  (see 
Section  5.4.2)  is  used  in  our  calculation.  In  a  two  plane,  par¬ 
allel  electrodes  design,  with  the  consumption  of  the  anode, 
the  cell  gap  will  increase,  and  thus  the  ohmic  resistance  will 
increase.  Larger  cell  gap  means  larger  ohmic  resistance.  At 
lower  cell  voltage,  the  current  density  is  higher,  and  the  cell 
gap  has  a  greater  effect  on  the  current  density. 

7.4.  Effect  of  parasitic  current  on  selectivity-corrosion 
inhibitor 

One  of  the  objectives  in  the  optimization  of  the  cell  per¬ 
formance  is  to  increase  the  selectivity  of  the  main  current 


Fig.  7.  Effect  of  cell  gap  on  current  density  at  different  constant  voltages  and  a  electrolyte  flow  rate  of  8  cm/s. 
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Fig.  8.  Effect  of  corrosion  current  on  selectivity  at  a  flow  rate  of  8  cm/s  and  cell  gap  of  0.2  cm  using  sodium  stannate  as  a  corrosion  inhibitor. 


density  over  parasitic  current  density.  From  Fig.  8,  one  can 
see  that  if  we  can  decrease  parasitic  current  by  a  factor  of  10, 
the  selectivity  will  be  much  higher.  In  using  pure  aluminum, 
and  no  corrosion  inhibitor  addition  in  the  electrolyte,  the 
selectivity  will  be  very  low  at  low  current  density  (main 
reaction)  (curve  2).  In  the  case  of  using  aluminum  alloy 
(A1  1199)  and  the  additive  sodium  stannate  (Na2SnC>3),  a 
corrosion  inhibitor,  in  the  electrolyte,  the  parasitic  current 
will  decrease  by  a  factor  of  ten  as  we  have  shown  in  our 
experimental  studies. 

7.5.  Effect  of  electrolyte  conductivity  on  current  density 

Fig.  9  shows  the  effect  of  electrolyte  conductivity  on 
current  density  at  different  constant  cell  voltages.  Small 
changes  in  electrolyte  conductivity  have  significant  effect 
on  the  current  density.  With  an  increase  in  conductivity,  the 


current  density  will  significantly  increase  in  all  three  cases. 
So  in  the  choice  of  the  electrolyte,  the  conductivity  will  be 
an  important  factor.  The  specific  electrolyte,  its  concentra¬ 
tion  and  temperature  all  affect  the  conductivity  of  the  elec¬ 
trolyte.  A  KOH  aqueous  solution  has  higher  conductivity 
than  a  NaOH  aqueous  solution  when  their  concentrations 
are  equal.  That  is  why  many  people  choose  KOH  aqueous 
as  the  electrolyte.  Of  course,  other  factors  should  also  be 
taken  into  account  such  as  the  effect  of  different  electrolyte 
on  the  crystallization  of  the  products  and  the  regeneration 
of  aluminum  by  the  industrial  Hall-Herout  process. 

7.(5.  Activation  overpotential  and  ohmic  loss  at 
different  cell  gaps 

Activation  overpotential  and  ohmic  loss  at  different  cell 
gaps  are  listed  in  Table  4.  The  concentration  overpotential 


Fig.  9.  Effect  of  conductivity  on  current  density  at  different  constant  cell  voltages  at  a  flow  rate  of  8  cm/s  and  cell  gap  of  0.2  cm. 
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Table  4 


Cell  Overpotential  and  ohmic  losses  at  0.8  V  and  8  cm/s  velocity 


Cell  gap  (cm) 

Current  density  (mA/cm2) 

Activation  loss  (V) 

Anode 

Cathode 

Ohmic  loss  (V) 

Total  loss  (V) 

0.2 

528 

1.32  (68.4%) 

0.47  (24.4%) 

0.14  (7.2%) 

1.93  (100%) 

0.4 

437 

1.25  (64.7%) 

0.45  (23.7%) 

0.23  (11.6%) 

1.93  (100%) 

0.6 

381 

1.20  (62.2%) 

0.43  (22.2%) 

0.30  (15.6%) 

1.93  (100%) 

0.8 

338 

1.16  (60.1%) 

0.42  (21.8%%) 

0.35  (18.1%) 

1.93  (100%) 

is  not  included  in  this  table  for  it  is  very  small  compared 
with  the  activation  overpotential  and  ohmic  loss.  At  a  flow 
rate  of  8  cm/s,  the  effect  of  gas  bubbles  on  ohmic  resistance 
is  small,  the  ohmic  resistance  is  considered  to  be  constant 
along  the  cell  channel.  When  the  cell  gap  increases  from 
0.2  to  0.8  cm,  the  anode  overpotential  decreasing  from  68.4 
to  60.1%  (8.3%  net  decrease);  the  cathode  overpotential  is 
relatively  stable,  it  decreases  from  24.4  to  21.8%  (2.6% 
net  decrease);  the  ohmic  loss  increases  from  7.2  to  18.1% 
(10.9%  net  decrease).  The  current  density  decreases  from 
528  to  338mA/cm2  (190mA/cm2  or  36.0%  net  decrease). 
Although  at  a  fixed  cell  voltage  of  0.8  V,  the  ohmic  loss  only 
increases  7.2-18.1%,  the  effect  on  current  density  is  signif¬ 
icant,  36%  current  density  will  be  lost  due  to  the  increase 
of  ohmic  loss  (or  cell  gap). 

Table  5  shows  overpotential  and  ohmic  losses  at  two  cell 
gaps,  different  cell  voltages.  We  can  see  that  the  ohmic 
loss  become  important  at  higher  current  density  (lead  to  the 
decrease  of  current  density  at  the  same  cell  voltage),  and 
selectivity  is  higher  at  higher  current  density. 

7. 7.  Gas  fraction  variation  with  cell  height 

Fig.  10  indicates  the  gas  fraction  variation  with  height  in 
the  aluminum-air  cell.  Calculations  represented  in  Fig.  10 
have  been  made  when  there  is  no  additive  in  the  electrolyte 


to  suppress  the  corrosion  reaction.  It  is  clear  that  the  volume 
fraction  takes  up  a  significant  percentage  of  the  channel  at 
the  exit  rising  from  20%  for  Case  1  to  as  high  as  40%  for 
Case  3.  These  volume  fractions  cause  large  reductions  in  the 
cell  conductivity  and  indicate  the  importance  of  additives 
under  these  conditions.  Fig.  11  clearly  depicts  the  impor¬ 
tance  of  increasing  the  electrolyte  velocity  to  sweep  out  the 
bubbles  generated  from  the  corrosion  reaction.  The  higher 
velocities  markedly  reduce  the  gas  volume  fraction  yielding 
higher  electrolyte  conductivity. 

If  we  used  a  corrosion  inhibitor  in  the  calculations  for 
Fig.  10  to  reduce  the  exchange  current  density  of  the  cor¬ 
rosion  reaction  to  /po  =  1.1  mA/cm 2  the  volume  fraction  at 
a  cell  height  of  17  cm,  flow  velocity  of  1  cm/s,  1.3  V  and 
200  mA/cm2  reduces  to  0.02  from  about  0.2.  The  additive 
causes  a  reduction  by  a  factor  of  10  in  the  volume  fraction. 
Other  parameters  are  unchanged  in  the  corrosion  kinetic 
equation. 

7.8.  Effect  of  gas  fraction  on  electrolyte  conductivity 

Electrical  conductivity  calculations  using  the  values  of 
Fig.  10  demonstrate  the  reduction  in  electrical  conducti¬ 
vity  in  the  bulk  electrolyte  flow.  The  reference  value  for 
the  electrolyte  conductivity  is  0.8  S/cm.  As  the  gas  frac¬ 
tion  builds  up  with  height  in  the  three  cases  illustrated, 


Table  5 

Overpotential  and  ohmic  losses  at  different  cell  gaps,  different  cell  voltages  and  a  flow  rate  of  8  cm/s. 
Cell  gap  —  0.2  cm 


Voltage,  V 

0.4 

0.6 

0.8 

1 

1.2 

1.4 

932.4 

715.9366 

529.05 

375 

254.53 

166 

5  m  /p 

85.85 

53.2843 

31.06 

16.8942 

8.5507 

4.0446 

Anode  loss 

1.5289 

1.4314 

1.3206 

1.195 

1.054 

0.8999 

Cathode  loss 

0.5566 

0.5098 

0.4685 

0.4335 

0.4048 

0.3822 

Ohmic  loss 

0.2405 

0.1848 

0.1369 

0.0975 

0.0672 

0.0439 

Total  loss 

2.326 

2.126 

1.926 

1.726 

1.526 

1.326 

Cell  gap  —  0.4  cm 

Voltage,  V 

0.4 

0.6 

0.8 

1 

1.2 

1.4 

732.9783 

577.7783 

440.075 

322.37 

226.408 

152.39 

5m /p 

55.5849 

36.3326 

22.4186 

12.9513 

6.964 

3.483 

Anode  loss 

1.4402 

1.3529 

1.2535 

1.1404 

1.0122 

0.869 

Cathode  loss 

0.5136 

0.4794 

0.4485 

0.4212 

0.3978 

0.3785 

Ohmic  loss 

0.3722 

0.2937 

0.224 

0.1644 

0.116 

0.0785 

Total  loss 

2.326 

2.126 

1.926 

1.726 

1.526 

1.326 

Bulk  electrolyte  gas  fraction,  f(y) 
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Fig.  10.  The  gas  fraction  variation  with  height  in  the  aluminum-air  cell  at  a  flow  rate  of  1  cm/s. 


Fig.  11.  The  gas  fraction  variation  with  height  in  the  aluminum-air  cell 
at  an  electrolyte  velocity  of  8  cm/s. 


the  electrical  conductivity  is  significantly  reduced  (Figs.  12 
and  13). 

7.9.  Anode  surface  concentrations 

Table  6  shows  the  surface  concentrations  on  the  anode  at  a 
cell  height  of  17  cm  and  /po  =  11  mA/cm2.  The  surface  con¬ 
centrations  of  OH-  and  Al(OH)4_  are  similar  at  electrolyte 
velocities  of  1  and  8  cm/s.  At  current  densities  between  520 
and  530  mA/cm2,  the  surface  concentration  of  OH-  is  2% 
less  than  the  concentration  in  the  bulk  electrolyte  and  the 
surface  concentration  of  Al(OH)4_  is  8%  larger  than  the 
concentration  in  the  bulk  electrolyte. 

Table  7  shows  the  surface  concentrations  on  the  anode 
at  a  cell  height  of  17  cm  and  zpo  =  1.1  mA/cm2.  The  sur¬ 
face  concentrations  of  OH-  and  Al(OH)4_  are  also  similar 
at  electrolyte  velocity  of  1  and  8  cm/s.  At  current  density 


Fig.  12.  The  reduction  in  electrical  conductivity  in  the  bulk  electrolyte  flow  at  electrolyte  velocity  of  1  cm/s. 
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Fig.  13.  The  reduction  in  electrical  conductivity  in  the  bulk  electrolyte  flow  at  electrolyte  velocity  of  8  cm/s. 


Table  6 

The  surface  concentrations  on  the  anode  at  cell  height  of  17  cm  and  ip o  =  llmA/cm2 


v  —  1  cm/s 

Cell  voltage,  V 

0.4 

0.6 

0.8 

1 

1.2 

1.4 

920 

705 

521 

368 

250 

163 

5m /p 

84.0154 

52.0133 

30.2559 

16.42 

8.2969 

3.9179 

COH,0/COH,b 

0.9606 

0.9727 

0.9815 

0.9879 

0.99 

0.995 

C  AlOH ,  0  /  C  AlOH ,  b 

1.1893 

1.1296 

1.0859 

1.0543 

1.033 

1.019 

v  =  8  cm/s 

Cell  voltage,  V 

0.4 

0.6 

0.8 

1 

1.2 

1.4 

932 

716 

529 

374.9092 

254.03 

165.9951 

5m/p 

85.849 

53.2843 

31.06 

16.8942 

8.5507 

4.0446 

COH,0/COH,b 

0.9606 

0.9723 

0.9812 

0.9876 

0.99 

0.995 

C  AlOH ,  0  /  C  AlOH ,  b 

1.1893 

1.1317 

1.0878 

1.0559 

1.033 

1.0197 

Table  7 

Surface  concentration  at  cell  height  of  17  cm  and  ip o 

=  1.1  mA/cm2 

V  =  1  cm/s 

Cell  voltage,  V 

0.4 

0.6 

0.8 

1 

1.2 

1.4 

911 

703 

522 

371.8 

253.2285 

165.53 

5m/p 

905.41 

548.13 

315.137 

170.083 

85.7438 

40.4854 

COH,0/cOH,b 

0.897 

0.9318 

0.9559 

0.9725 

0.9834 

0.9903 

CAlOH,0/CAlOH,b 

1.502 

1.3324 

1.2143 

1.133 

1.0797 

1.0457 

V  =  8  cm/s 

Cell  voltage,  V 

0.4 

0.6 

0.8 

1 

1.2 

1.4 

917 

707.55 

524.556 

372.86 

253.7788 

165.844 

5m /p 

894.07 

545.58 

315.029 

170.33 

85.9372 

40.5864 

COH,0/COH,b 

0.911 

0.9393 

0.9599 

0.9745 

0.9844 

0.9908 

C  AlOH ,  0  /  C  AlOH ,  b 

1.4326 

1.2951 

1.1948 

1.1236 

1.075 

1.0435 

520-530  mA/cm2,  the  surface  concentration  of  OH-  is  4% 
less  than  the  concentration  in  the  bulk  electrolyte  and  the 
surface  concentration  of  Al(OH)4-  is  20%  larger  than  the 
concentration  in  the  bulk  electrolyte.  Comparing  the  results 
of  Tables  6  and  7,  we  can  see  that,  at  high  corrosion  rate, 
the  surface  concentration  polarization  can  be  relieved. 


8.  Conclusions 

We  have  modeled  the  aluminum-air  cell  performance 
using  the  equations  we  have  developed  in  this  study.  The 
model  prediction  of  cell  performance  shows  good  agree¬ 
ment  with  experimental  data.  For  better  cell  performance, 
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our  model  studies  suggest  the  use  of  higher  electrolyte  flow 
rates,  smaller  cell  gaps,  higher  conductivities,  lower  par¬ 
asitic  current  densities,  and  operation  at  moderate  current 
densities.  From  our  analysis,  we  have  determined  for  rela¬ 
tively  high  flow  rates  high  (>8  cm/s)  in  an  aluminum  air  cell, 
the  gas  effect  on  current  density  distribution  is  small.  We 
can  also  conclude  that  only  the  activation  and  ohmic  over¬ 
potentials  are  important  in  the  model  calculations.  Further 
the  proper  additive  can  increase  the  selectivity  by  a  factor 
of  10. 
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